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Summary  and  status  overview 


This  effort  is  focused  on  fundamental  investigations  of  mixing,  chemical-reaction, 
and  combustion  processes,  in  turbulent,  subsonic,  and  supersonic  free-shear  flows. 
The  program  is  comprised  of  an  experimental  effort;  an  analytical,  modeling,  and 
computational  effort;  and  a  diagnostics  and  data- acquisition-development  effort. 
The  computational  studies  are  focused  on  fundamental  issues  pertaining  to  the 
hydrocarbon-ignition/-combustion  and  numerical  simulation  of  compressible  flows 
with  strong  fronts,  in  both  chemically-reacting  and  nonreacting  flows. 


1.  Accomplishments  and  new  findings 


To  investigate  the  effects  of  various  parameters  on  turbulent-shear-layer  mixing 
and  combustion  at  high  Reynolds  number,  experiments  were  conducted  to  ex¬ 
plore  the  influence  of  inflow  conditions.  Chemically-reacting,  incompressible,  shear- 
layer  flows  were  investigated  in  the  GALCIT  Supersonic  Shear  Layer  Facility,  at 
freestream  conditions  U2/U1  ~  0.4  and  P2/P1  —  1  and  a  local  Reynolds  number. 
Res  =  pAU6(x)/n  ~  2  x  10®. 


Fig.  1  Normalized  temperature-rise  data.  Left:  untripped  boundary  layers.  Right: 
tripped  high-speed  boundary  layer.  Reactant  compositions  -  diamonds:  ^  = 
8;  triangles:  (f>  =  1/8;  asterisks:  (j)  =  1/8,  reduced  chemical-kinetic  rate. 


3 


Chemically-reacting  “flip”  experiments  allow  us  to  deduce  the  structure  and  amount 
of  molecular-scale  mixing,  information  unobtainable  from  direct  scalar  imaging  at 
such  high  Reynolds  numbers.  Two  sets  of  temperature-rise  data  from  such  ex¬ 
periments  are  presented  in  Fig.  1.  The  only  difference  between  these  is  a  0.8  mm- 
diameter  trip  wire  on  the  splitter  plate  (high-speed  side),  50  mm  upstream  of  its 
trailing  edge.  The  data  are  measured  far  downstream,  at  xfQx  ~  3300,  where  x 
is  the  streamwise  coordinate  and  6i  is  the  high-speed  boundary  layer  momentum 
thickness,  and  at  a  large  value  of  the  pairing  parameter,^  P  ~  47.  By  these  criteria, 
the  flow  may  be  regarded  as  fully-developed.^’^  The  change  in  inflow  conditions  can 
be  seen  to  have  a  signiflcant  effect  on  all  scales  of  the  flow:  the  shear-layer  growth 
rate,  8/x^  has  decreased  by  21%;  the  mixed-fluid  fraction  (efficiency  of  molecular- 
scale  mixing),  has  increased  by  11%;  and  the  mixed-fluid  composition  ratio  has 
decreased  by  9%.  Additionally,  the  data  indicate  a  change  from  a  non-marching 
scalar  probability- density  function  (pdf),  to  a  marching  pdf,  when  the  boundary 
layer  was  tripped.  These  observations  suggest  a  shear-layer  behavior  that  depends 
not  only  on  local-flow  properties,  but  also  on  upstream  conditions,  reminiscent  of 
results  from  low-dimensionality,  nonlinear  (chaotic)  systems.'^ 


Fig.  2  Transverse  jet  at  Rej  ~  1.0  x  10®  and  Uj/Uoo  —  34.  Left:  LIF  image  in  plane 
parallel  to  free-stream  (log-intensity  scale).  Right:  Slice  perpendicular  to 
free-stream. 

Experiments  on  the  development  and  far-field  structure  of  a  liquid-phase  turbulent 
jet  issuing  into  a  cross-flowing  uniform  stream  were  conducted  at  a  jet  Reynolds 
number,  Rej  =  Ujdj/i/  ~  1.0  x  10®,  for  a  velocity  ratio,  Uj/Uoo  —  34.  Laser-induced 
fluorescence  (LIF)  images  measured  scalar-species  concentration  in  two  orthogonal 
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Fig.  3  Shape  complexity  pdf  of  islands  and  lakes.  Solid  line:  power-law  fit  (in  these 
coordinates). 

slices  of  the  transverse  jet  (Fig.  2).  The  data  show  wake  structures  similar  to  those 
reported  previously,^  e.g.^  “fingers”  of  dyed  jet-fluid  extending  from  the  tunnel  wall 
to  the  main  body  of  the  jet.  These  are  transported  downstream  with  nearly  uniform 
velocity  after  being  formed  near  the  jet  exit. 


As  part  of  our  investigation  of  the  geometry  of  scalar  isosurfaces  in  liquid-phase 
turbulent  jets,  an  analysis  of  area-perimeter  properties  of  isoscalar  “islands”  and 
“lakes”  derived  from  2-D  images  of  jet-fluid  concentration  images  has  been  con¬ 
ducted.  The  images  are  slices  in  a  plane  normal  to  the  jet  axis  in  the  far  field  and 
correspond  to  a  fluid  Schmidt  number  of  Sc  ~  2.0  x  10^  and  flow  Reynolds  num¬ 
bers  of  Re  4.5  X  10^,  9.0  x  10^,  and  18  x  10®. For  such  data,  scalar  level  sets 
form  disjoint  “islands”  and  “lakes”,  depending  on  whether  the  interior  is  at  a  lower, 
or  higher,  scalar  level,  respectively.  Island /lake  statistics,  such  as  size  and  shape 
complexity,  are  important  in  many  contexts:  for  chemical  reactions  and  combustion 
in  nonpremixed  hydrocarbon  turbulent  flames,  for  example,  combustion  is  largely 
confined  to  the  instantaneous  stoichiometric  (isoscalar)  surface.®  In  2-D,  such  a 
measure,  0,2  ■,  dubbed  shape  complexity ,  can  be  defined  as. 


1  <  O2 


2(nA)'/^ 


<  00  , 


with  P  the  perimeter  and  A  an  island/lahe  area,  and  (f^2)iniii  =  ^  attained  for  a  cir¬ 
cle.  Corresponding  extensions  can  also  be  made  for  higher-dimensional  embedding 


spaces. 
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Analysis  of  the  liquid-phase  jet  data  described  above  indicates  that  a  power  law  for 
over  3  decades  in  size  (6  decades  in  area),  or,  equivalently,  log-Poisson  statistics, 
provides  a  good  approximation  for  the  pdf  of  shape  complexity,  with  similar  be¬ 
havior  expected  in  turbulent  flows  in  general.  Such  distributions  are  necessary  in 
modeling  burning-time  distributions  in  nonpremixed  combustion,  for  example.® 

Numerical-simulation  investigations  were  conducted  on  the  ignition  characteris¬ 
tics  of  hydrocarbon  fuel  blends,  at  conditions  relevant  to  high-Mach-number,  air- 
breathing-propulsion  vehicles.  A  two-point-continuation  method  was  employed, 
with  a  detailed  description  of  molecular  transport  and  chemical  kinetics,  focus¬ 
ing  on  the  effects  of  fuel  composition,  reactant  temperature,  additives,  and  imposed 
strain  rate.  It  captured  the  entire  S-curve  that  describes  the  processes  of  vigorous 
burning,  extinction,  and  ignition.  The  results  demonstrate  that  ignition  of  such 
fuel  blends  is  dominated  by  the  synergistic  behavior  of  CH4  and  C2H4.  A  fuel  tem¬ 
perature  Tfuei  =  950  K  was  used  in  the  simulations.  For  low  air  temperatures,  e.g., 
Tair  1050  K,  addition  of  small  amounts  of  CH4  enhances  C2H4  ignition,  while  for 
Tair  >  1200  K,  CH4  inhibits  ignition.  The  ignition  of  blends  of  CH4  and  C2H4  was 
significantly  promoted  through  independent  additions  of  small  amounts  of  H2,  and 
F2  with  NO.  The  latter  yields  F-radicals  which  effectively  consume  hydrocarbon 
molecules.  An  important  finding  of  these  continuing  investigations  is  that  the  GRI 
2.1  Cl,  C2  kinetics  mechanism,  which  was  developed  to  describe  CH4  combustion, 
cannot  represent  ethylene  (C2H4)  ignition  and  kinetics.  Interestingly,  it  is  missing 
a  key  species,  absent  which,  the  strain  rate  sufficient  to  prevent  ignition  is  underes¬ 
timated  by  a  factor  of  25-30.  A  first  report  of  this  work  was  presented  at  the  recent 
International  Combustion  Symposium  (Boulder,  CO).^° 

Developments  on  Image  Correlation  Velocimetry  (ICV),  a  method  that  extracts 
velocity-field  information  from  sequences  of  images  of  a  convected  scalar,^ ^^2  ^re 
continuing.  This  has  been  used  to  analyze  two-dimensional  slices  of  low-speed  liquid- 
phase  flows,  with  both  particles  and  fluorescent  dyes,  as  well  as  spacecraft  imaging 
of  the  Jovian  atmosphere.  An  important  extension  to  a  three-dimensional  imple¬ 
mentation  is  presently  in  progress,  with  envisaged  initial  testing  on  existing  three- 
dimensional  laboratory  data  on  turbulent  jets  {e.g.,  Ref.  13),  as  well  as  a  variety 
of  direct  numerical-simulation  (DNS)  data.  An  application  to  Ozone  Tagging  Ve¬ 
locimetry  image  data,^'*  is  also  in  progress,  in  collaboration  with  R.  W.  Pitz  and 
his  group  at  Vanderbilt  University,  who  are  developing  this  technology.  Finally,  in 
collaboration  with  C.  M.  Vagelopoulos  and  F.  N.  Egolfopoulos  (USC),  an  extension 
of  the  original  particle-streak  velocimetry  technique,^®  targeted  for  strained  flames 
and  hydrocarbon  laminar  flame-speed  measurement  applications. 
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mixing  at  high  Reynolds  numbers:  effects  of  inflow  conditions,”  J.  Fluid  Mech. 
(accepted).  GALCIT  Report  FM98-1. 
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Slessor,  M.  D.,  Zhuang,  M.  &  Dimotakis,  P.  E.  1998  “Turbulent  shear-layer 
mixing:  growth-rate  compressibility  scaling,”  GALCIT  Report  FM98-9. 


5.  Interactions/Transitions 

Visits/interactions/participation/presentations  at  meetings,  conferences,  seminars: 

-  Bond  (presenter),  C.  L.,  “Mixing  in  high  Reynolds  number  shear  layers.”  GAL¬ 
CIT  Fluid  Mechanics  Seminar  (6  March  1998). 

-  Catrakis  (presenter),  H.  J.,  and  Bond,  C.  L.,  “Scale  distributions  of  geometric 
transects  of  level  sets,”  contributed  presentation.  Am.  Phys.  Soc.,  Div.  Fluid 
Dynamics,  annual  meeting  (San  Francisco,  CA,  November  1997). 

-  Dimotakis  (presenter),  P.  E.,  “Turbulence,  fractals,  and  mixing.”  Plenary  lec¬ 
ture,  Am.  Phys.  Soc.,  Div.  Fluid  Dynamics,  annual  meeting  (San  Francisco, 
CA,  November  1997). 

-  Papalexandris  (presenter),  M.  V.,  Leonard,  A.,  and  Dimotakis,  P.  E.,  “Unsplit 
Schemes  for  Multi-Dimensional  Systems  of  Hyperbolic  Conservation  Laws  with 
Source  Terms,”  contributed  presentation.  Am.  Phys.  Soc.,  Div.  Fluid  Dynam¬ 
ics,  annual  meeting  (San  Francisco,  CA,  November  1997). 

-  Shan  (presenter),  J.  W.,  Laidlaw,  D.  H.,  Gornowicz,  G.  G.,  Lang,  D.  B.,  and  Di¬ 
motakis,  P.  E.,  “Three-Dimensional  Space- Time  Structure  of  Turbulent  Jets,” 
contributed  presentation.  Am.  Phys.  Soc.,  Div.  Fluid  Dynamics,  annual  meet¬ 
ing  (San  Francisco,  CA,  November  1997). 

-  Slessor  (presenter),  M.  D.,  Bond,  C.  L.,  and  Dimotakis,  P.  E.,  “The  effect  of  ini¬ 
tial  conditions  on  turbulent  shear  layers,”  contributed  presentation.  Am.  Phys. 
Soc.,  Div.  Fluid  Dynamics,  annual  meeting  (San  Francisco,  CA,  November 
1997). 

-  Dimotakis  (presenter),  P.  E.,  “Image  correlation  velocimetry  and  the  flow  over 
an  accelerating  NACA-0012  airfoil.”  GALCIT  Fluid  Mechanics  Seminar,  with 
R.  Henderson  (9  January  1998). 

-  Slessor  (presenter),  M.  D.,  “Dynamics  and  Mixing  of  Turbulent  Shear  layers.” 
GALCIT  Fluid  Mechanics  Seminar  (6  March  1998). 


11 


-  Dimotakis  (presenter),  P.  E.,  Catrakis,  H.  J.,  Cook,  A.  W.,  and  Patton,  J. 
M.,  “On  the  geometry  of  two-dimensional  slices  of  irregular  level  sets  in  tur¬ 
bulent  flows.”  Invited  lecture,  2”**  Monte  Veriia  Colloquium  on  Fundamental 
Problematic  Issues  in  Turbulence  (22-28  March  1998,  Ascona,  Switzerland). 

-  Dimotakis  (presenter),  P.  E.,  and  Slessor,  M.  D.,  “Mixing  in  high  Reynolds 
number  shear  layers.”  Plenary  lecture,  AIAA  Fluid  Dynamics  meeting  (Albu¬ 
querque,  NM,  June  1998). 

-  Egolfopoulos  (presenter),  F.  N.,  and  Dimotakis,  P.  E.,  “Effects  of  additives 
on  hydrocarbon  ignition  at  high  strain  rates.”  International  Combustion 
Symposium  presentation/paper  (Boulder,  CO,  August  1998). 

-  Clavin  (U.  Marseilles),  Paul;  Discussions  on  detonation,  supersonic  test  facili¬ 
ties,  and  related  research  issues  (August  1998  visit  to  Caltech). 


Consultative  and  advisory  functions  to  other  labs,  agencies,  especially  AF  and  DoD 
(include  institutions,  locations,  dates,  names): 

-  Dimotakis,  P.  E.,  Lawrence  Livermore  National  Laboratories.  Consulting  on 
compressible  turbulence,  inertial-confinement  fusion,  high-fluence  laser-doubling 
crystal  growth  (1993  to  present). 

-  Dimotakis,  P.  E.,  Member  of  NAS/NRC  Committee  that  reviewed  DOE’s  In¬ 
ertial  Confinement  Fusion  (ICF)  program  and  the  National  Ignition  Facility 
(work  completed  August  1997). 

-  Dimotakis,  P.  E.,  Member  Mitre/ JASON  group. 


Transitions  (cases  where  knowledge  resulting  from  this  research  is  used  or  will  be 
used  in  a  technology  application): 

-  None 


New  Discoveries,  inventions,  patents: 


None 


